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ABSTRACT 


A coherent bistatic vegetation scattering model, based on a 
Monte Carlo simulation, is being developed to simulate 
polarimetric bi-static reflectometry at WHF/UHF-bands 
(240-270 MHz). The model is aimed to assess the value of 
geostationary satellite signals of opportunity to enable 
estimation of the Earth’s biomass and root-zone soil 
moisture. An expression for bistatic scattering from a 
vegetation canopy is derived for the practical case of a 
ground-based/low altitude platforms with passive receivers 
overlooking vegetation. Using analytical wave theory in 
conjunction with distorted Born approximation (DBA), the 
transmit and receive antennas effects (i.e., polarization, 
orientation, height, etc.) are explicitly accounted for. Both 
the coherency nature of the model (joint phase and 
amplitude information) and the explicit account of system 
parameters (antenna, altitude, polarization, etc) enable one 
to perform various beamforming techniques to evaluate 
realistic deployment configurations. In this paper, several 
test scenarios will be presented and the results will be 
evaluated for feasibility for future biomass and root-zone 
soil moisture application using geostationary 
communication satellite signals of opportunity at low 
frequencies. 


Index Terms— Bistatic, UHF, root-zone soil moisture, 
biomass, signal of opportunity, discrete scatterer. 


1. INTRODUCTION 
In the past decade, there has been increasing interest in 


using signals of opportunity (SoOp) for remote sensing of 
Earth [1]. Most of these past studies have focused on 


collection or modeling of Global Navigation Satellite 
Systems (GNSS) signal reflections over ocean surface for 
estimating wind direction and speed or altimetry. On the 
other hand, application of GNSS-Reflectometry technique 
to land is rather scarce and relies mostly on experimental 
observations [2-4]. Sensitivity to soil and vegetation by 
GNSS’s bistatic signals at L-band has been, to some extent, 
experimentally demonstrated and have been also supported 
by few theoretical studies [5-9]. However, L-band signals 
are limited to penetrate through dense vegetation and deep 
into the soil column. Thereby, GNSS-based approaches can 
provide estimate of surface soil moisture up to a depth of 
~0.05 meter through moderate amounts of vegetation. In 
recent past, to address this issue, a new SoOp methodology 
is suggested, where the signals transmitted by the Military 
Satellite Communication (MilSatCom) System’s UHF (240- 
270MHz) Follow-On program) are utilized to measure 
properties of reflecting targets by recording reflected signals 
using a simple passive microwave receiver [10, 11]. 

The low frequency SoOP approach has the potential of 
measuring root-zone soil moisture (top 1-m) through heavy 
vegetation. These measured quantities would provide 
information critical for understanding global water 
resources and the vertical moisture gradient responsible for 
evapotranspiration and runoff. In this study, a low 
frequency (VHF/UHF) bi-static electromagnetic scattering 
model is being developed to supplement recent 
experimental efforts and to exploit further their potential. 
Since the low-frequency SoOP experiments are so far 
carried out from either ground based or airborne platforms, 
we will restrict the modeling study to such geometries, but it 
can be fairly easily extended to the receivers in low Earth 
orbit. The model will be used to conduct sensitivity 
analysis and to determine optimum system configuration for 


estimating root-zone soil moisture and biomass parameters 
for the Earth surface covered with heavy forest/vegetation. 

The antenna characteristics and orientation play a key 
role in received signal. For ground-based systems, the 
antenna radiation pattern projected on the surface is not 
uniformly distributed in phase, amplitude, or polarization. 
The bistatic radar coefficient calculated by the model thus 
needs to have the antenna characteristics embedded in it as 
well as the statistical and physical properties of the terrain to 
mimic the real measurement setting. Previous studies have 
usually assumed plane wave illumination/scattering [5-7], 
and ignored such antenna effects. In the following section, 
a brief description of the model will be provided for 
calculating the response of the antenna due to the coherent 
bistatic scattered field from a vegetated surface. 


2. FORMULATION 


In this paper, a general bistatic scattering case is considered, 
where two antennas do not have their main beam axis 
pointing at each other and both of which are overlooking the 
vegetation. One antenna system is associated with the 
signal of opportunity transmitter while the other refers to the 
passive receiver. The vegetation is represented as ensemble 
of canonical scatterers located above a flat ground. The 
transmitter is located far away. The local incident angles 
are assumed to be constant (parallel) and spreading loss 
effects due to the incoming wave are ignored, but the 
spreading loss and sphericity of the scattered wave are 
considered due to proximity of the platforms that operate 
close to the ground. This configuration will represent 
“plane wave incidence” and “spherical wave scattering”. 

It is important to realize that antennas cannot be 
constructed to produce pure polarization states. There will 
always be some nonzero cross polarization level that could 
degrade the measurement. Additionally, the polarization 
state of the antenna is usually defined along the main beam 
direction. The polarization, however, varies for increased 
off-axis angles over the radiation pattern. Polarization 
states of the antennas are taken into account through a 
polarization (field) rotation matrices. As a result, the 
formulation considers variations of both the strength and 
polarization states of the received wave along the beam 
direction. 

Applying ‘discrete scatter’ modeling approach [12], the 
forest layer is replaced by a slab of equivalent homogenous 
medium by the Foldy theory [13]. An expression for the 
scattered field due to a single particle immersed in the slab 
of mean medium over a flat ground surface is derived. The 
received signal is then expressed as a weighted coherent 
sum involving the antenna spatial spectrum and bistatic 
scattering. The results of the single particle scattering are 
summed over all types and particles by taking into account 
of the polarization mixing (rotation matrices), the antenna 
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Figure 1: 4-channel VHF receiver. 


impurities (antenna pattern matrices with co- and cross-pol 
components), and the spreading loss. 

The received field has three main terms: (1) direct term, 
(2) specular term, (3) diffuse term. The first one represents 
the line of sight or the shortest path between the antennas 
while the second term denotes the forward scattering along 
the specular direction, where the incoming and _ the 
scattering waves make the same angle with respect to the 
surface normal. Lastly, the diffusely scattered waves arrive 
at the receiver antenna from a wide range of angles in both 
azimuth and elevation due to scattering from the illuminated 
volume. An average diffuse term is obtained by a sufficient 
number of realizations of vegetation through Monte Carlo 
simulations. 

Each term is a “complex” field quantity and thus 
includes both amplitude and phase information. This 
enables one to decompose the field into a orthogonal 
polarization basis in order to simulate polarimetric bi-static 
reflectometry for various polarization combinations. 


3. SIMULATION SETTINGS 


In the simulation, several scenarios will be tested for the 
recently developed 4-channel VHF receiver hardware [11]. 
The instrument is shown in anechoic chamber for 
measurement of radiation characteristics in Fig. 1. Each 
channel receives VHF signal picked up by the linear dipole 
tuned to frequency range 240-270 MHz. Due to broad 
beam dipole patterns, each channel will receive the direct as 
well as reflected signal (reflected from the vegetation) 
transmitted by the MilSatCom satellites. The instrument is 
capable of recording complex field quantities that allows to 
apply beamforming. Using digital beam forming in the post 
processing, the direct signal and the reflected signals are 
separated for soil moisture and root zone soil moisture 
extraction. 


Zenith, direct signal, horizon, and reflected signal 
directions will be considered for antenna orientations. For 
each configuration, the 4 channel array beam will be steered 
into direct and reflected direction to isolate the signals. The 
model will be simulated for considering various 
combination of system and target parameters. Among the 
antenna parameters, effects of antenna height, antenna 
cross-polarization, and beamwidth width will be analyzed to 
understand platform and antenna requirements for such 
applications. For each system configuration, the target 
parameters such as soil moisture, vegetation type, height, 
and biomass level, will be varied to conduct sensitivity 
analysis. 


4. CONCLUSION 


Knowledge of the Earth biomass, root zone soil moisture 
(including surface soil moisture) are important earth science 
parameters because of their strong influence on the Earth’s 
climate system. Currently, these Earth science parameters 
are measured using airborne and space borne microwave 
radiometers/scatterometers operating at L-band (1.4 GHz) 
frequencies. However, the L-band frequency range is 
inadequate for measuring these parameters for the Earth 
surface covered with heavy forest/vegetation. In the recent 
past, NASA has shown great interest in using lower 
frequency band such as VHF/UHF for these parameter 
measurements. In this study, newly developed coherent 
bistatic vegetation scattering model, based on a Monte Carlo 
simulation, will be utilized to simulate polarimetric bi-static 
reflectometry at VHF/UHF-band. Preliminary performance 
of the model and sensitivity analysis by the model will be 
presented. 
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